Splitting of the Dipole and Spin-Dipole Resonances 



Sam M. Austin 1 *, E. Adamides 1 , A. Galonsky 1 , T. Nees 1 , W.A. Sterrenburg 1 ^ , 
D.E. Bainum 2 *, J. Rapaport 3 , E. Sugarbaker 4 ^, C.C. Foster 5 , CD. Goodman 5 , 
D.J. Horen 6 **, C.A. Goulding 7tt , and M.B. Greenfield 7 " 
1 National Superconducting Cyclotron Laboratory and Department of Physics and Astronomy 
Michigan State University, East Lansing, MI 48824 
2 Department of Physics, Emporia State University, Emporia, KS 66801 
3 Department of Physics, Ohio University, Athens, OH 4 5701 
4 Department of Physics, University of Colorado, Boulder, CO 80309 
5 Indiana University Cyclotron Facility, Indiana University, Bloomington, IN 47405 
6 Oak Ridge National Laboratory, Oak Ridge, TN 37830 
7 Department of Physics, Florida A&M University, Tallahassee, FL32307 

(February 8, 2008) 

Abstract 

Cross sections for the 90 ' 92 ' 94 Zr(p, n) reactions were measured at energies of 
79.2 and f 19.4 MeV. A phenomenological model was developed to describe the 
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variation with bombarding energy of the position of the L = 1 peak observed 

in these and other (p,n) reactions. The model yields the splitting between 

the giant dipole and giant spin dipole resonances. Values of these splittings 

are obtained for isotopes of Zr and Sn and for 208 Pb. 
PACS numbers: 25.40.Kv, 24.30.Cz 
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I. INTRODUCTION 



Inelastic scattering and charge exchange reactions induced by intense neutrino fluxes 
during supernova explosions play an important role in nucleosynthesis in supernovae. For 
example, they may produce rare light isotopes such as 7 Li, n B and 19 F |IJ that are not 
sufficiently produced by the usual processes of nucleosynthesis and may modify the distri- 
bution of elements formed during the r-process [[J. Observation of supernova neutrinos in 
terrestrial detectors sensitive to neutrino flavor || may allow one to determine or limit the 
masses of the /x and r neutrinos and to determine whether oscillations occur. From the point 
of view of nuclear physics, all these phenomena share a common ingredient: they depend on 
the L = allowed (Gamow Teller) and L = 1 forbidden strengths of transitions induced by 
the various neutrino flavors. It may be possible to measure some of these strengths directly 
using neutrino beams, but it will be necessary to infer most of them from the cross sections 
for hadronic charge exchange and inelastic scattering reactions. 

Unfortunately, information on the nature of the L = 1 strength is not yet sufficiently 
detailed for such applications. The spectrum observed at low momentum transfer in a typical 
(p, n) reaction is dominated by three isovector excitations: the isobaric analog state (IAS), 
the giant Gamow Teller Resonance (GTR), and the giant dipole (L = 1) resonance || [|. 
It is well known that the IAS contains almost all the Fermi strength and that the Gamow 
Teller peak contains about 60% of the sum rule value. Much less is known about the dipole 
peak and it exhibits a rather unusual behavior; it moves to lower excitation energy as the 
bombarding energy increases as shown in Fig. 1. 

The qualitative reason for this behavior seems clear; the L = 1 peak contains two separate 
structures: the giant dipole resonance (GDR) formed by a spin transfer of S — 0, and the 
spin dipole resonance (SDR) formed by S — 1. These resonances lie at different positions, 
but are sufficiently wide that they overlap. Because their relative excitation is expected 
to change in the energy range from 40 to 200 MeV, that of the SDR becoming relatively 
stronger at higher energy 0, the L = 1 peak will shift toward the SDR as the energy 
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increases. The amount of the shift depends on the energy splitting of the GDR and SDR; 
it should then be possible to obtain an estimate of this splitting from an analysis of the 
observed shifts. 

The Zr isotopes are an excellent place to study the shift phenomenon. There have been 
0, n) studies on 90 Zr §-0] at 41, 45, 120, and 200 MeV, and on other Zr isotopes at 45 MeV 
|| and 200 MeV |lTJ. This paper describes measurements of the 90 ' 92 ' 94 Zr(p, n) reactions at 
two intermediate energies, 79.2 and 119.4 MeV; a model that can account for the energy 
variation observed for the most complete data set ( 90 Zr); and the application of this model 
to extract the splitting of the S = and S = 1 components of the L = 1 resonance for 
several other nuclei. The basic assumptions of this model are similar to those employed by 



Taddeucci et al. |12| in the description of L = charge exchange reactions with S = and 
5=1. 

The model also provides an improved estimate of the position of the GDR obtained from 
(p, n) reaction studies. In earlier work it has simply been assumed that at bombarding 
energies near 50 MeV the observed excitation is dominated by the GDR; we show that this 
is not the case. We note that one cannot simply obtain the position of the GDR in the 
charge exchange channel from photonuclear studies on the target nucleus. In nuclei with 
substantial values of N-Z, the dominant excitation in (p, n) is the T — 1 component of the 
GDR where T is the isospin of the target nucleus. The collective shift of this resonance 
differs from those of the T and T + 1 components seen in the target. Phenomenologically, 



this is accounted for by an isotensor component of the effective isospin potential HJT3|] , which 



is poorly known. Consequently, the location of the T — 1 GDR must be taken as unknown. 

II. EXPERIMENTAL PROCEDURES AND RESULTS 

Data were taken using the beam swinger time of flight facility at the Indiana University 
Cyclotron Facility |TJj], with two 15 cm x 15 cm x 100 cm detectors, located 46.5 m and 



70.15 m from the target. Bombarding energies were 79.2 MeV and 119.4 MeV, nominally 
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80 and 120 MeV. Targets were enriched 90 > 92 > 94 Zr metal with thicknesses of 77.5, 25.4, and 
41.6 mg/cm 2 , respectively. The overall resolution was between 415 and 540 keV (580 and 
675 keV) at 79.2 MeV (119.4 MeV), significantly smaller than the width of the L = 1 peaks 
we are concerned with here. Spectra were taken at about 0.1, 4, 8, 12, 16 and 19.5 deg(lab). 
Those for 90 Zr at 4 deg are shown in Fig. 2; at this angle the IAS and the GT resonances, 
which peak at deg, and the L = 1 resonances, which peak at 7-8 deg are visible. 

Spectra were fitted with sums of Gaussian peaks plus a quadratic background, as shown 
in Fig. 2. The resulting fits are satisfactory. For angles of 4 deg and greater a single gaussian 
provides a good description of the L = 1 strength. This is somewhat surprising because the 



S = 0, S = 1 splitting is expected to be in the range of 3.0 to 4.5 MeV [p~5|| , not so much 
smaller than the width of the GDR, and because the relative contributions of the GDR and 
SDR are expected to change from 1:1 to about 1:20 over the energy range studied. 

A numerical experiment was performed to evaluate the uncertainties that might arise 
from this procedure. Gaussian peaks with widths of 4.5 MeV (representing the GDR) and 
8 MeV (representing the SDR), and 3.9 MeV apart, were added and their sum was fitted 
by a single Gaussian. This was done for peak ratios from 1:1 to 1:20. In no case did the 
centroid of the fitted Gaussian differ from the actual centroid by more that 0.15 MeV, and 
most differences were much smaller. In another test, the centroid of the counts above the 
fitted background was calculated and compared to the centroid of the fitted Gaussian. The 
difference was 0.17 MeV. This gives one some confidence in the procedure followed at about 
the ±0.1 MeV level. The resulting energies for the dipole peak are shown in Fig. 1; the 
L = 1 peak shifts to lower energies by about 2 MeV as the bombarding energy changes from 
40 to 200 MeV. The peak positions are consistent with centroid energies obtained by others 
at both lower and higher bombarding energies (see Fig. 4), lending some confidence to the 
procedures used here. 
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III. MODEL AND RESULTS 



We outline here a phenomenological model that describes the observed shift with bom- 
barding energy in terms of the positions of the GDR and the SDR. The basic idea is as 
outlined in the Introduction: that the energy of the L = 1 peak moves toward the energy of 
the SDR as the bombarding energy increases. The notation we use is shown in Fig. 3. We 
denote the energy of the 5 = component, the J w = 1~ GDR, by Eq. Similarly, the energy 
of the 5 = 1 component, the SDR consisting of 0~, 1~, and 2~ states has a combined cen- 
troid denoted by E\. The states with different J will have different distributions of strength 
with excitation energy. We lack the information to treat these differences in excitation in 
more detail, but believe that the effects are probably small compared to other uncertainties 
in the model. 

Because of the nature of the two-nucleon interaction mediating the transitions, the 
strength of the 5 = excitation of the GDR is comparable to that of the 5 = 1 excitation 
of the SDR at 45 MeV, but 5=1 becomes much stronger than 5 = as the bombarding 
energy increases toward 200 MeV. One then expects the observed centroid of the L = 1 
strength to move toward E\ with increasing bombarding energy. Rather general arguments 



1~5| indicate that the 5 = GDR lies at higher excitation energy than the centroid of the 



5=1 SDR so this motion corresponds to a decreasing excitation energy. Detailed the- 
oretical calculations [0] qualitatively confirm this expectation. See Refs. f§-|6| for general 
discussion and references. 

The position, C, of the centroid of the L = 1 excitations (including both the GDR and 
SDR) at a bombarding energy E p is given by 

c= °oEo + °iE 1 =Eo __* 1 /*o A 
o- + o"i 1 + 0-1/0-0 

where A=Eq - E\ and o"o(o"i) is the cross section for 5 = 0(5 = 1) transfer. One can express 
C in terms of Eq and A if one knows the ratio o"i/o"o- We assume here that for the L = 1 
excitations this ratio is directly proportional to the ratio of the forbidden 5 = 1 and 5 = 
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beta decay strengths, as is observed for the L = case ]12[ where the analogous excitations 



are the Gamow Teller and IAS resonances. This assumption is confirmed for L = 1 by the 
DWIA calculations of Gaarde et al. ||16|| . Then for small momentum transfers one has [12j 

Jn 



K 



£?(SDR) 

fi(GDR) ' 1 ' 



Jr 

where J aT and J T are the volume integrals of the parts of the effective interaction medi- 
ating S = 1 and S = transitions respectively, the -B's are the corresponding transition 
strengths and the K's are kinematic factors whose ratio is very near 1.0 for our case. Np 
= (N aT /N T )L = i where the iV's are distortion factors which differ from each other by about 
15% at the first peak in the L = 1 angular distribution. For L = transitions N^\J aT / J T \ 2 
= a 2 E 2 where = (N aT / 'N T ) L=0 , a = (55 MeV) -1 and E p is the bombarding energy |12 



This relationship is accurate over the energy range of 45 to 200 MeV studied here. One then 
obtains 

aa 2 El , s 

C = E - A . 3 

1 + aa 2 E* 

Here a is given by (A^/iV^)[i?(SDR)/I?(GDR)]. In principle, experimental results can 
be used to fix a if data at three energies are available (there are three parameters A, E 
and a). Fitting the 90 Zr results yields a value of a near 1, but with a large uncertainty due 
to the limited accuracy of the centroid data. It is then desirable to obtain an estimate of a 
at least partly from other considerations. 

The ratio of the B's depends on both the sum rule strengths for the transitions and 
the fraction / of those strengths (quenching factor=/) found in the low excitation giant 
resonance peaks of interest here. Estimates with uncorrelated and with RPA ground states 
I]J give 5(SDR)/5(GDR) = 3.25 and 3.2, respectively. Gaarde, et al. [|TIJ have also shown 



that the effect of RPA correlations is small. For the 5=1 quenching factor we take / ~ 
0.5, comparable to that obtained experimentally for the L = GT resonance Q. A value 
of / in the range of 0.3 to 0.5 has been observed by Gaarde et al. |TjJ for masses of 40 to 



208. Several authors, e.g. Drozdz et al. ||18| , have shown that perhaps 30% of the total 



L = 1 strength is moved to high energies. Strength at higher excitation has also been seen 
in polarization transfer experiments and in multipole decompositions of (p, n) and (n, p) 
data Our results then apply to the resonant part of the cross section, not that part 

lying higher in the continuum. Implicitly the model assumes that the largest concentration 
of strength, the peak in the excitation energy distribution for L = 1, gives a reasonable 
estimate of the position of the collective strength for comparison with models. There is 



some indication in Ref. |18| that the position of maximum strength does not change greatly 
when continuum effects are considered. For S = the quenching factor is presumed to be 
near 1.0, in analogy with the L = 0, 5* = (Fermi transition) case. 

Finally N^/N^ is obtained from the DWIA calculations evaluated at the peaks of the 
relevant angular distributions (e.g. near 4.5° for L = 1 and at 0° for L = at 200 MeV on 
Pb) and typically has a value near 1.0. Based on all these considerations we take a = 1.5 
± 0.5 as covering the values of interest. 

Fitting Eq. 3 with a — 1.5 ± 0.5 to the data for 90 Zr yields the results shown in Fig. 4. 
Similar procedures were followed for other isotopes for which data are available at widely 
separated bombarding energies. For 92 ' 94 Zr the excitation energy data are shown in Fig. 1; 
for 112 - 11 6,i24g n data are available a t 45 and 200 MeV [|,|TT]]; and for 208 Pb data are available 



at 41,45,120 and 200 MeV The results are shown in Table I. Also shown in Table I 



is the position of the GDR measured in pion induced charge exchange reactions. These data 
are not in good agreement with the (p, n) results for 90 Zr. However, the difference between 
the results of single and double charge exchange reactions for 208 Pb indicates that there may 
be systematic uncertainties not yet accounted for. 

Krmpotic, Nakayama and Pio Galeao |15 and Auerbach and Klein have discussed 



the nature of the giant first forbidden resonances and have given values for the energies of 
the GDR and SDR: their results are given in Table I. They are in qualitative agreement 



with the results of a detailed calculation by Bertsch, Cha, and Toki |21]] which places the 
centroid of the SDR near 18 MeV, below the GDR. Theory and experiment agree closely for 
208 Pb and within the uncertainties for the Zr isotopes, although the theoretical results are 
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smaller. The agreement is less good for the Sn isotopes, perhaps because experimental data 
are available only at 45 and 200 MeV. It is difficult to estimate the theoretical uncertainties, 
but they could be a few MeV, mainly in the position of the spin-dipole excitation ||22| . 



IV. DISCUSSION AND SUMMARY 

It is perhaps surprising that the the experimental and theoretical values of Eq and A 
shown in Table I are in reasonably good agreement, since the experimental results do not 
measure the position of the high-lying continuum strength and the theoretical estimates 
include all strength. We do not have an explanation for this result, except possibly that the 
position of maximum strength calculated in RPA does not change much when continuum 
effects are included; see Ref. [fTB| . Calculations of the positions of the resonance peaks would 



be more appropriately compared with the present experiments. 

Information similar to that in Table I can in principle be obtained from a comparison of 
the positions of the SDR observed in high energy (p, n) reactions and the GDR observed in 
pion charge exchange; however, the energy shift approach outlined here at the least provides 
an independent measure of the GDR-SDR splitting, and arguably is less subject to system- 
atic uncertainties arising from different experimental approaches. The present procedure is 
based on simple assumptions and simple Gaussian shapes for the dipole resonances. The 
procedure could be greatly improved if data of sufficient statistics and angular range were 
available at several energies for each target so that subtraction or multipole fitting proce- 
dures could be used to isolate L = 1 strength. With better data, it would become worthwhile 
to model corrections owing to the (relatively small) contributions of the higher isospin (T) 
components of the excitation; here we have implicitly neglected all but the T = Tq — 1 
components. 

Even in its present form the model provides a quantitative explanation of the shift of 
the L = 1 excitation with bombarding energy, a measure of splitting of the spin-dipole and 
dipole resonances, and an improved estimate of the position of the GDR peak (see Table I 
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and Fig. 4). The results appear in a general way to support the assumptions made by Fuller 
et al. H in their calculations of the the response of a flavor-dependent supernova neutrino 
detector. They assume that the spin-dipole resonance does not lie at higher excitation energy 
than the GDR, in agreement with our finding. However, the present results only apply to 
the part of the strength that is in the resonant peaks, while the response of the detector will 
depend also on the part of the strength that lies in the continuum at higher E x . In addition, 
it will be necessary to use a model to relate the observed splitting in the charge exchange 
channel to that in 208 Pb. 

The results also indicate that past estimates of tensor contributions to the isospin split- 
ting of the giant dipole resonance (e.g. Sterrenburg et al. ||) must be modified. These 
estimates assumed that the GDR was at the position of the L = 1 peak observed for bom- 
barding energies near 50 MeV. Our results show that it lies at higher E x , reducing the values 
of E_ = Et — Et-i obtained in that paper. 

To summarize, cross sections for the 90 > 92 > 94 Zr(p, n) reaction have been measured at en- 
ergies of 79.2 and 119.4 MeV. The positions of the observed L = 1 excitation have been 
extracted and compared with those at other energies; the excitation strength is observed to 
move toward lower excitation energy as the bombarding energy increases. A simple phe- 
nomenological picture was developed to describe the observed energy variation and was used 
to extract the splitting of the giant dipole (GDR) and giant spin-dipole (SDR) resonances. 
The results are in agreement with theoretical models as to the scale of the splitting, but the 
value of A we obtain is usually somewhat larger than predicted. 

We have had valuable conversations with N. Auerbach, G. Bertsch, D.M. Brink, 
K. Nakayama and H. Toki. E. Adamides acknowledges the financial support of the Academy 
of Athens-Greece from the V. Notara Bequest. This research was supported by the National 
Science Foundation and the Department of Energy. 
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FIGURES 

FIG. 1. Energies of the L = 1 resonance for the Zr isotopes measured from the IAS; values 
are from the present work, and from Refs. [§-ll|. The uncertainties of the present data include 



an estimate of the systematic errors; this also appears to be the case for the data from the other 
references. 

FIG. 2. Spectra from 90 Zr(p, n) at 79.2 and 119.4 MeV, with background used in the analysis. 
The energy is measured from the IAS. 

FIG. 3. Diagram showing assumed model and definition of terms. The zero of the energy scale 
is the energy of the IAS. 

FIG. 4. Comparison of the L = 1 resonance positions for 90 Zr and Eq. 3. Energies are measured 
from the IAS. The abscissa is the bombarding energy and the curves are non-linear least square 
fits to the data. 



13 



TABLES 

TABLE I. Energy a of the Dipole (E ) and of the Dipole-Spin-Dipole Splitting (A) 





Present 


b 


Theory 




(vr+,vr°),(7r+,7r-) 


Target 


Eq 


A 


Eq 


A 


Eq 




MeV 


MeV 


MeV 


MeV 


MeV 


90 Zr 


16.7±1.2(0.6) 


4.0±1.3(0.9) 


15.3, 14.6 


3.0 


14.2±0.5 d 


92 Zr 


17.1±1.6(0.8) 


5.1±1.7(1.1) 


15.0 


3.4 




94 Zr 


16.9±1.9(1.0) 


5.4L2. 1(1.4) 


14.7 


3.8 




112 Sn 


16.6±2. 1(1.0) 


6.8±2.7(1.7) 


13.8 


2.8 




116 Sn 


17.5±2.5(1.2) 


8.6±3.1(1.8) 


13.1 


3.3 




124 Sn 


14.4±2. 2(1.0) 


7.5±2.8(1.7) 


12.3 


4.3 




208p b 


11.0±1.5(0.8) 


4.7±2.0(1.4) 


10.3, 9.9 


4.7 


10.2±0.9 C , 8.6±0.5 f 



a Energies are with respect to the isobaric analog state 

b Calculated for a = 1.5. The uncertainties include an allowance for the uncertainty of ± 0.5 in a: 
one half the difference in £o and A for fits with a = 1.0 and 2.0 is added linearly to the fit error. 
The first number following the ± sign is the total error; the fit error is given in parentheses. 



From Ref. [15|, except the second values of Eq for 9U Zr and Pb are from Ref. [p0| 



d Mean of results in Refs. ||§|25| 
e Ref. |H (7r+,7r-) 
f Rcf. ||, (vr+,7r°) 
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